Dipolar bosons on an optical lattice ring 
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We consider an ultra-small system of polarized bosons on an optical lattice with a ring topology 
interacting via long range dipole-dipole interactions. Dipoles polarized perpendicular to the plane 
of the ring reveal sharp transitions between different density wave phases. As the strength of the 
dipolar interactions is varied the behavior of the transitions is first-order like. For dipoles polarized 
in the plane of the ring the transitions between possible phases show pronounced sensitivity to the 
lattice depth. The abundance of possible configurations may be useful for quantum information 
applications. 

PACS numbers: 67.85.Hj, 03.75.Lm 



The behavior of atoms in optical lattice potentials has 
attracted a lot of attention in the past decade [l|. The 
impressive degree of control over atomic samples gives 
unprecedented chances for studying models coming from 
condensed matter physics, e.g. the Hubbard model for 
fermions or the Bose-Hubbard model for bosons J2J]- This 
allows an investigation of relatively large systems which 
are needed for quantum phase transition problems [3j 
in order to minimize finite size effects. A subsequent 
progress reaches the stage of single atom or single site 
imaging [J| . Also optical lattices of practically arbitrary 
configurations can be prepared |5|,[6| opening up new pos- 
sibilities. 

While ultra-cold atom-atom interactions are typically 
dominated by contact terms, long range dipole-dipole in- 
teractions are becoming of great interest - for recent re- 
views see [3,la|. Experimental realizations of condensates 
of Cr [9( and very recently of dysprosium [lOj open up 
a pathway for novel physics of atomic species exhibiting 
strong dipolar interactions. On the other hand, owing to 
tunable Feshbach resonances one is able to control scat- 
tering lengths between particles very precisely [ll| ■ This 
tuning may help also to enhance the role of dipolar in- 
teractions. 

Quite impressive work has been performed on polar- 
ized dipolar gases in one-dimensional lattice systems (see 
e.g. [12h17I |). Here by using the density matrix renormal- 
ization group (DMRG), quantitative studies have been 
carried out efficiently even for relatively large systems. 
New phases such as density waves (absent for purely on- 
site interactions as in the Bose-Hubbard model) can then 
be observed [l^[l3J]. Longer range interactions allow for 
the presence of the Haldane insulating phase (with a non- 
trivial so-called string order [Ij, LLJI)- In these works it 
is assumed that the dipoles are polarized perpendicularly 
to the one-dimensional lattice so the dipolar interaction 
is repulsive. This way one can avoid possible instability 
due to "head-tail" attraction between dipoles. 

The importance of the dipoles' orientation has been 



nicely visualized in the study of a dipolar condensate 
placed in a toroidal trap \lj$. Assuming the dipoles to 
be oriented in the plane of the ring it becomes apparent 
that the dipole-dipole interactions become non-uniform 
along the ring, with repulsive interactions where locally 
the ring is perpendicular to the polarization axis and at- 
tractive where the ring reinforces "head-tail" alignment. 
Using the mean field Gross-Pitaevskii description it has 
been shown [18| that sufficiently strong dipole-dipole in- 
teractions split the condensate into two parts. Another 
study [19( considered very few bosons or fermions polar- 
ized at an arbitrary angle with respect to a small ring 
by using exact diagonalization. In both these works no 
additional optical lattice was assumed. 

The results presented in this paper combine the two 
situations. We consider a few bosons on an optical lat- 
tice realized with a ring topology and consider dipolar- 
interactions with different polarizations. Interestingly, 
this leads to qualitatively novel behavior of the model. 

Assuming the lattice to be sufficiently deep so that the 
excited "bands" may be neglect ed 3 (for limitations on 
that approximation see e.g. |2Cl. l2l|). the system may be 
described by an extended Bose-Hubbard model: 
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(1) 
The first part of the Hamiltonian is just the usual one di- 
mensional Bose-Hubbard model, where b\ and bi are the 
lattice creation and annihilation operators and hi = b]bi, 
is the corresponding number operator. The coefficients J 
and U correspond to the strength of hopping and contact 
on-site interaction terms. 

The last term in the Hamiltonian describes the dipole- 
dipole interactions which provide coupling between all 
the sites. The coefficient Vij depends both on the mutual 
orientation between two sites and the distance between 
them: 
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where lUj(r) is a Wannier wavefunction localized on site 
i, fio is the permeability of free space, and fj, is the mag- 
netic moment of the specific particle (to be specific we 
consider magnetic interactions, although a similar anal- 
ysis holds for electric dipoles). 8 is the angle between 
the axis of polarization and the vector connecting two 
dipoles r — r' where the dipoles are positioned at r and 
r'. Since the interactions are anisotropic, they can be 
either attractive or repulsive depending on the relative 
position of the dipoles. It is exactly here where the ring 
topology of the lattice makes a difference - otherwise, the 
Hamiltonian, Eq.Jl}, is identical to the one-dimensional 
models considered previously (la . KL4| . It is also worth 
stressing that we consider dipole interactions between all 
sites while large scale DMRG treatments were typically 
limited to the next-nearest neighbors couplings. 

The on-site i — j contribution of the dipole-dipole in- 
teraction is assumed to be incorporated into the U term, 
while for i ^ j we assume that the lattice is sufficiently 
deep that the angle 8 and distances may be measured 
with respect to the centers of the sites (i.e. taken out 
of the integrals). Defining the lattice constant a as the 
distance between nearest neighbors a = |R, — R*+i| we 
may construct dimensionless positions of sites along the 
ring R; = Rj/o. Then 
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where V measures the strength of the dipole-dipole in- 
teractions and 8ij corresponds to the angle between the 
orientation of the dipoles and the vector connecting them 
at sites i and j. 

Consider first the case of dipoles polarized in a direc- 
tion perpendicular to the plane of the ring. The angle 
is always 8ij = ir/2 and the dipoles repel each other and 
the magnitude of the interaction depends only on the 
distance between sites via the 1/r 3 law. In this situation 
the model becomes quite similar to the one-dimensional 
linear models with periodic boundary conditions (with 
quantitative differences coming from different Vij depen- 
dence on the distance i — j due to a different geometry) . 
Contrary, however, to typical DMRG studies of large lin- 
ear systems [131 fl4| we shall consider 8 bosons sitting on 
8 sites (i.e. with unit mean density). This allows us to 
use standard exact diagonalization routines [22j . 

Of course for such a small system we should not speak 
about quantum phase transitions but merely consider 
crossovers between different ground state arrangements 
(note that quite a small dipole system was discussed re- 
cently for a triple well potential [23). Surprisingly, stan- 
dard measures (such as a superfluid fraction or a con- 
densate fraction) used to describe large systems are also 
applicable and useful for our case. At the same time finite 
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Figure 1: (color online) Superfluid fraction f s (solid black 
lines), condensate fraction f c (red dashed lines) both rep- 
resented on left-side vertical axis as a function of V/U for 
U = 1,4,7,10 (panels (a) to (d), respectively). Right axis 
scale corresponds to green dot-dashed lines that represent 
variance of occupation on different sites. For small V/U we 
observe an increase of the superfluid (and condensate) frac- 
tion with V both in the SF phase (for U — 1 panel (a) as well 
as for U = 10 panel (d)) where the increase of V causes tran- 
sition from MI to SF phase. Around V/U = 0.5 a transition 
to a density wave state [2, 0] occurs, more sudden for larger 
U, as manifested by a rapid drop of f s and f c . For larger 
V/U sharp transitions between different density waves occur 
as manifested in f c as well as in variance of occupation A. 



size effects become important and interesting. The con- 
densate fraction, f c , is obtained from the largest eigen- 
value of the one-body density matrix (OBDM). Observe 
that for our system f c > 1/8. Following Roth and Bur- 
nett [24( the superfluid fraction, f s , is obtained by includ- 
ing Peierls phase factors in the Hamiltonian to simulate 
the slow motion in the system (in our case this corre- 
sponds to a spinning of the ring). 

The Hamiltonian, ([I]), depends on three parameters, 
J, U, and V. We assume J = 1 thus expressing U, V, 
and H in units of J. 

Fig. [T] presents / s , f c , as we ll as the variance of the 
occupation, A = -\/(rt, 2 ) — (n-i) 2 as a function of the ratio 
of V and U . Different panels refer to different values 
of U. In particular, panels (c) and (d) correspond to 
values of U at which the infinite system at V — would 
be deep in the gapped Mott insulating phase. A rough 
explanation of the behaviour of f s and f c at small V/U 
can be obtained assuming that the hopping amplitude is 
negligible, J « 0. At V — a gap AE = U separates 
the ground state energy from the first excited state. In 
this situation one would expect f s = f c = 0, and the 
observed non-zero values are due to finite-size effects. A 
simple calculation shows that for small values of V/U the 
gap is reduced to AE = U — V, which entails an increase 
in the superfluid and condensate fraction. 

This increase in the superfluid and condensate fraction 
is not observed in panel (a), corresponding to a value of 
U/J at which the infinite system would be superfluid 



in the absence of dipolar interactions. These quantities 
start off already at their maximum value, and decrease 
as V is increased. Panel (b) refers to U/J — 4, which 
would be slightly inside the Mott insulating state for an 
infinite lattice at unitary filling and in the absence of 
dipole interactions. Indeed the behaviour of the finite 
system under examination is more alike to that observed 
at larger U/J, in panels (c) and (d), although not as 
sharply peaked. 

In panels (b)-(d) the SF fraction exhibits a sudden de- 
crease with increasing dipole interaction at V/U ~ 0.5, 
and the system enters into a density-wave state charac- 
terized by a site occupation number alternating between 
and 2. We denote this situation as [2, 0] state. As 
we discuss in the following, different density-wave states 
take over at larger values of V/U. As can be seen from 
the plots in panels (b)-(d), the (sharp) crossover points 
between subsequent density-wave regimes are essentially 
independent of the value of U, since at negligible J/U 
the only relevant parameter is V/U. 

Since all sites are equivalent, the system is invariant 
under rotation by 2ir/M where M = 8 is the number of 
sites. The ground state shares this symmetry. There are 
two almost degenerate lowest lying states corresponding 
to even and odd combinations of [2,0] and [0,2] states 
satisfying (n/) — 1 at every site. To identify the den- 
sity wave state, we construct the density matrix of the 
ground state and trace out all sites except one. This 
yields the probability distribution P(n) of bosons at a 
single site (see Table H}. The second column (for which 
the V/U value corresponds to a [2, 0] density wave phase) 
shows that n = and n — 2 are occupied with the same 
probability being almost 1/2 confirming the alternating 
empty and doubly occupied sites, the remaining n values 
have non vanishing occupation due to finite tunneling. 
Observe that just below V/U = 3 (see fig. [3J>) a jump 
appears in f c and the variance A - a transition occurs 
to the density wave state in which two sites are occupied 
by 3 bosons and one by 2 bosons, as it is evident from 
P(n) distribution in the third column of Table fl] The 
next jump in A indicates a transition to one site occu- 
pied by 4 bosons and two by 2 bosons while the last jump 
creates two sites with 4 bosons with the rest remaining 
empty. Eventually for much higher V/U (not shown) 
one observes that all bosons accumulate at a single site, 
leaving 7 sites empty. This is because accumulation of 
particles at one site costs little energy (proportional to 
U - a vanishingly small parameter in this limit) while 
occupation of any two sides gives a nonzero term propor- 
tional to V. This behaviour is in a sharp contrast with a 
cristal-like structure obtained for a ring topology without 
a lattice [19J where the dipole interaction enhances repul- 
sion due to a contact term thus small distance between 
bosons are avoided. Recall, however, that we assume 
that the dipole-dipole interaction between particles at a 
single site is incorporated in U term already, so there is 
no additional energy costs for particles to accumulate at 
a single site. 
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Table I: Probability distributions at a given site for different 
ratio of dipolar to on-site coupling, V/U for [7 = 4 revealing 
different density wave arrangements. 



Naturally, the types of density waves observed are 
strongly system size dependent. For 9 bosons on 9 sites 
we observed a state with one site occupied by 3 bosons 
and three sites by 2 bosons, then a crossover into a repe- 
tition of [3, 0, 0] sequence and for large V/U one site with 
5, another with 4 atoms and the remaining sites empty. 
It should be noted that some density waves taking over 
the ground state of the system appear because truly long- 
range interactions are considered. If, as chosen in |l2l4l4| , 
the range of the interaction is truncated, some density- 
wave configurations never crop out. For instance, if no 
interactions are included beyond next-nearest-neighbors, 
the configuration where all bosons accumulate at a single 
site cannot appear, no matter how strong the interaction 
is. Importantly observe that crossovers between different 
density waves resemble phase transitions of the first kind 
(the jump in A and also in / c ). 

Consider now the the orientation of dipoles in the plane 
of the ring. In this case interactions between two sites 
are not only dependent on the distances separating them 
but are also influenced by their relative positions on the 
lattice. Due to the small size of our system (we again 
consider 8 bosons on 8 sites) different orientations of the 
dipole polarization with respect to the octagon are possi- 
ble, the results depend quantitatively on this orientation. 
Further on we choose the most symmetric case when the 
polarization axis coincides with two edges of the octagon 
(compare Fig. |2K)- Then, from symmetry considerations 
we have two types of sites: 4 located at the end of the 
edges that are parallel to the polarization axis (type-I) 
and 4 located at the end of the edges perpendicular to 
that axis (type-II). Due to the symmetry, sites of a given 
type are equivalent. Recall that the dipole-dipole inter- 
action terms in the last sum in Eq. ([1]) can be either 
repulsive or attractive. As we will discuss in the follow- 
ing, the physics of the system is dominated by the strong 
attractive interaction between dipoles occupying neigh- 
bouring type-I sites. 

Fig. [2}d shows system properties for a shallow lattice 
U / J = 1. The increase of V/U has the effect of transfer- 
ring particles from sites of type II to sites of type I. For 
V/U ps 0.5 a change of the system's properties appears. 
The SF fraction drops from values close to unity to al- 
most zero; the variance of the particle number in type-I 
sites, A/, exhibits a sudden and significant increase; sim- 
ilarly a drop in the CF, f c , is observed, while a second 
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Figure 2: (color online) A ring shaped lattice (a), the arrow in- 
dicates polarization of dipoles. Due to symmetry sites 2,3,6,7 
are equivalent (type-I), similarly 1,4,5,8 (type-II). Panel (b) 
shows results for a shallow lattice, U/J — 1. The solid black 
line shows a SF fraction with a drop around V/U — 0.5. This 
drop is accompanied by similarly rapid increase of the num- 
ber variance A/ (blue dash-dotted line, right axis scale) as 
well as changes of condensate fraction, f c (dashed red line) 
and the second eigenvalue of OBDM, g c (dotted red line). On 
the other hand the mean occupation n/ of type-I sites (blue 
dashed line, right axis scale) increases slowly and smoothly in 
the whole studied range of V/U. For the discussion see text. 



eigenvalue of the OBDM becomes large. This suggests 
that, as V/U exceeds 0.5 the system essentially consists 
of the superposition of two coherent halves, each localized 
at a pair of neighbouring sites of type I. This superpo- 
sition of "macroscopic"states is typical of a "double well" 
system where an effective attractive interaction is domi- 
nant [23, [26j , and gives rise to the sudden increase in the 
number fluctuations observed in Fig. [2Jd. 

This phe nomenology is reminiscent of the observations 
of Ref. [18| . where a situation similar to ours, but with- 
out the optical lattice, is analyzed adopting a Gross- 
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Table II: Probability distributions at sites of type I and II on 
either side of the V/U = 0.5 transition for U/J = 1, revealing 
ground state character change. 



Pitaevskii mean-field description. Similar to our case, an 
increase in the dipolar interaction causes at first a trans- 
fer of population around "hour 3 and hour 9", i.e. regions 
where dipoles arranged head to tail attract each other. 
As the interaction is further increased a symmetry break- 
ing - typical of attractive nonlinear systems - occurs, and 
the entire boson population localizes either around "hour 
3" or "hour 9". This symmetry broken situation is the 
mean-field counterpart of the macroscopic superposition 
occurring at the quantum level. The above description 
is confirmed by the inspection of the one site occupation 
probability distribution, P(n). Table II presents P(n) 
for values of V/U slightly below and slightly above 0.5. 
For type-II sites P{n) changes smoothly. On the other 
hand P(n) for type-I reveals that the for V/U > 0.5 the 
state becomes a superposition of the vacuum with a state 
centered around larger n. 

For larger values of the dipolar interaction (V/U > 30) 
the state of the system essentially consists of the su- 
perposition of two Fock states, -^|0,4, 4,0, 0,0, 0,0) + 

-4=|0,0, 0,0, 0,4,4, 0), whose condensate fraction is sig- 
nificantly smaller than 1/2. 

The "double-well" behavior becomes even more spec- 
tacular for intermediate values of U/J as demonstrated in 
Fig. [31 A sharp crossover is observed around V/U ~ 0.4 
for all the quantities plotted in panels (a)-(c). Like in 
Fig. [2] this separates the situation where a single coher- 
ent state is present from the situation where the system 
consists of the superposition of two coherent halves. The 
sharpness of the crossover may be linked to an extremely 
narrow avoided level crossing involving the ground-state 
of the system. Further increase of V/U smoothly depop- 
ulates type-II sites while the mean population of type-I 
sites becomes 2. At even larger V/U the ground-state 
eventually turns into the superposition of the two Fock 
states discussed above (as manifested by the near degen- 
eracy of the two largest eigenvalues of the OBDM, i.e. 
fc ~ 9c in Fig. [3]). We observe also that, at small values 
of V/U, the population transfer from type-II to type-I 
sites is accompanied by an increase in the superfluid and 
condensate fractions, as it is evident from panels (a) and 
(c). This effect has the same origin as explained when 
discussing Fig. [TJ At the transition point, apart from 
sharp changes in superfluid and condensate fractions, we 
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Figure 3: (color online) Intermediate U/J behavior. SF frac- 
tion is shown in panel (a) for U/J = 4 (black dotted line), 
U/J — 6 (red dashed line) and U/J — 8 (solid green line). 
Panel (b) presents the mean occupation and the variance of 
type-I states for U/J = 6 while panel (c) shows the V/U 
dependence (again for U/J = 6) of the condensate fraction, 
f c , i.e. the largest eigenvalue of the OBDM, and the second 
largest eigenvalue, denoted by g c . The sudden drop of SF 
fraction occurring at V/U = 0.4 in panel (a) is accompanied 
by a sharp increase of the occupation variance in type-I sites 
A/ [red dashed line in (b)]. The occupation ni (black line) 
shows there a small kink too. At the same V/U = 0.4 the 
sharp drop in ft is observed, while gb shows an increase. 



observe numerically also a small kink in the occupation 
of type-I sites (with, of course, a similar, reversed kink in 
occupation of type-II states) pointing to the redistribu- 
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Figure 4: (color online) Occupation of type-I sites, ni for 
U/J = 10 (solid black line) and U/J = 20 (solid red line). 
The corresponding variance Aj is plotted for U/J — 10 (black 
dashed line) and U/J = 20 (red dashed line). For a discussion 
see text. 



tion of densities between type I and type two sites due to 
a competition of tunneling and interactions. This effect 
is no longer present for stronger interactions (compare 
Fig.lU). 

For still stronger on-site repulsion U/J > 10, Fig. 21 
the system remains in an insulator-like state with van- 
ishing SF fraction for all V/U values. In this regime 
the state of the system is controlled by the competi- 
tion between the on-site repulsive and long-range dipo- 
lar interaction, and essentially consists of the super- 
position of a few Fock states. With reference to the 
U/J = 20 case in Fig. 5J the first plateau in the 
plotted quantities corresponds to the usual Mott state, 
1,1,1,1,1,1,1,1). Around V/U = 0.3 this is replaced 
by ^|0,1, 1,0,1, 2,2,1) + ^|1, 2,2,1, 0,1, 1,0). Next, 

around V/U = 0.4 the state i|0, 0,0, 0, 1, 3, 3, 1) + 
-4=|1, 3,3,1,0,0,0,0} takes over, which is then replaced 
by ^|0,0,0,0,0,4,4,0) + ^|0,4,4,0,0,0,0,0) around 

V/U = 0.65. We observe that the long-range character 
of the dipolar interaction is not essential for the above- 
described behaviour. It is the anisotropic character of 
dipolar interactions with attractive interactions between 
bosons in type-I sites which is responsible for the ob- 
served behaviour. Therefore the same qualitative feature 
would be observed also when a truncation in the range 
of interactions is included. 

To summarize we have shown that a simple model of 
a few bosons on a ring-shaped optical lattice can show 
very rich phenomena. At low V/U, dipole-dipole interac- 
tions tend to enhance the SF character. For dipoles po- 
larized perpendicularly to the ring plane various density 
wave states have been observed with sharp transitions 
between them as V/U was varied. For dipoles polarized 
in the plane of the ring the separation into two disjoint 



parts has been observed (as for the toroidal trap without 
the lattice [lM|). The presence of the optical lattices cre- 
ates additional sharp transitions due to rearrangement 
of dipoles in the parts of the ring where they attract 
each other. As J may be varied changing the depth of 
the lattice and the contact interactions strength may be 
modified using Feshbach resonances, the arrangements 
discussed should be accessible experimentally. 

This opens up a possibility for constructing quantum 
computer schemes in analogy to that realized for trapped 
ions (271 . [28fl . The latter, due to a strong Coulomb re- 
pulsion, are arranged with one ion per trap. The sys- 



tem studied may allow for controlled interactions between 
atoms by moving them between the sites in a controlled 
manner. 
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